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Abstract. In this article we address the problem of Euler's buckling instability 
in a charged semi-flexible polymer that is under the action of a compressive 
force. We consider this instability as a phase transition and investigate the role of 
thermal fluctuations in the buckling critical force. By performing molecular dynamic 
simulations, we show that the critical force decreases when the temperature increases. 
Repulsive electrostatic interaction in the finite temperature is in competition with 
thermal fluctuations to increase the buckling threshold. 
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1. Introduction 

Microtubules and actin filaments are important biopolymers that construct the 
mechanics of the eukaryotic cells [H [21 |3]. Microtubules the most rigid biopolymers 
of the cellular cytoplasm, organize to form asters and mitotic spindle inside the cells. 
The dynamics of these structures are essential in most of the physical processes of the 
cell cycle [11 15]. The actin filaments, the other less rigid biopolymers of the cell with 
rigidity of about 0.01 times that of the microtubules, are also play various roles in the 
cell [6]. The functioning of these filaments in the cell includes, the stability of the 
cell shape, crawling and the cell movements. Understanding the overall mechanical 
structure of a cell in terms of the behavior of the individual constituent polymers is of 
prime importance in cell mechanics. In addition to the bending properties, the buckling 
instability of the biopolymers in response to the pushing forces is being also observed 
[H IE]. In the case of microtubules two types of the forces initiate the buckling, the 
forces due to the polymerization when an astral microtubule reaches the cell cortex and 
the forces associated with the active molecular motors those are binding the different 
microtubules to each other or binding an individual microtubule to the cell cortex [QlfTO]. 
Actin filaments also experience the buckling instability because of the interaction with 
miosin motors |6]. 
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Figure 1. Schematic of a polymer that is under the action of a compression force /, 
at one of its ends. To investigate the stability of this polymer under this compressional 
force, we impose a transverse force /x, at the midpoint of the polymer. 



The buckling of the biopolymers are subject to either equihbrium or non equihbrium 
fluctuations. Equilibrium fluctuations due to the thermal equilibration with the ambient 
cytoplasmic fluid and the non equilibrium processes of the binding and unbinding of 
molecular motors to the filaments, enhance the buckling transition. In addition to 
the fluctuations, the other important effect that can change the buckling point, is the 
electrostatic interactions in a charged biopolymer. 

These observations show that a complete understanding of the mechanics at the 
cellular scale needs a detail knowledge of the mechanical properties of an individual 
filament. There are some efforts in considering the buckling of either charged or neutral 
filaments taking into account the fluctuation effects [HI Il2l [131 El IlSl HH]. Here by 
using the ESPRESSO package [17], we perfom the molecular dynamic simulations and 
focus on the numerical simulation for the buckling instability of either a charged or 
neutral filament that is subject to equilibrium thermally fluctuations. 

The rest of this article is organized as follows: In section 2 we introduce the model 
with its free energy and physical parameters. In section 3 some known analytic mean 
field results for neutral and charged polymers are summarized. We present the numerical 
results for a fluctuating filament in the section 4. 

2. The model 

Let us consider a stiff polymer with contour length L and bending modulus k embedded 
in a solution (Fig. [1]). The space configuration of this polymer is characterized by the 
position vector r(s) of a general points on the chains with contour distance s from one 
end. There is a uniform distribution of electric charges on this polymer. The free energy 
of this system has two important counterpart, the bending energy which is due to the 
chemical bonds between the adjacent monomers and the electrostatic repulsion between 
charges. The free energy of this system reads: 

^ = f r ■ £ + ^ r r '^'^^'^^^^'^'^ - ^^^'^^^ + 

where G*^^* is the contribution to the free energy due to the external forces acting on 
the chain and V{y{s) — r(s')), represents the electrostatic repulsion between the charges 
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on the polymer. The electrostatic interaction is screened by the mobile counter ions in 
the solution and is given by [T8] : 



2 p-r/\D 

= (2) 

where q represents the charge of the ions on the chain that are distributed uniformly 
with a minimum distance lo, e stands for the dielectric constant of the solution and Xd 
is the Debye-Huckel screening length. The dynamics of this chain in the over damped, 
low Reynold's number condition is given by the following equation: 

(3) 

where Q is the free energy density and the details of the hydrodynamics interactions 
in a simplified model are collected in two friction coefficients and C,±- Here is the 
friction coefficient for a small and cylindrical segment of the chain moving along the 
local tangent i, and denotes the corresponding friction for the motion along n, the 
unit vector that is locally normal to the polymer. The effects of thermal fluctuations 
are modeled by a random force ^{s,t) with Gaussian distribution. This means that the 
correlation function for this noisy force reads: 

{as,t)as',t'))=2D6is-s')6it-t'), (4) 

where D is the diffusion constant for a small segment of the chain. 

The aim of this article is to investigate numerically the response of a chain to a 
compressive force acting to one end of the chain. Before considering this problem we 
first present the results of some mean field descriptions. 



3. Mean field results 



We apply a constant compressive force f = — /x to one end of the polymer. This includes 
the following amount of mechanical energy to the free energy: 

G<=-* = -f .r(s = L), (5) 

minimizing this free energy allows us to study the equilibrium shape of the chain. This 
can be done for a set of different boundary conditions. Let us denote the angle that the 
local tangent vector at contour length s makes with x— axis by 6{s). 

Two important boundary are the following cases. Case (I): both ends of the 
chain are clamped (^(0) = 6{L) = 0) and case (II): both ends can freely rotate 
(0(0) = 9{L) = 0). Here the dot symbol represents the derivative with respect to contour 
length s. Physically the free rotating end is the case where no external torque is applied 
to the chain. Considering an electrically neutral chain in the mean field approximation 
where the effects of thermal fluctuations are not important, one can easily show that 
the equilibrium shape of the chain is the solution to the following differential equation: 

^eis) + fsme = 0. (6) 
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The straight chain solution with 9{s) = is a trivial solution that can exist for small 
forces. But one can see that for forces greater than a critical force /°, the above solution 
is not stable and a second solution with buckled shape for the chain initiates. The 
buckling threshold depends on the specific choice of the boundaries. For the condition 
where the two ends are clamped, case (I), the bucking force reads: /° = and for 
case (II) the buckling instability occurs at /° = 

Taking into account the effects of thermal fluctuation of a neutral semi-fiexible 
chain, enters a new length scale Lp = K/ksT, that is the persistence length of a 
fluctuating semi-flexible polymer. Saddle point analysis for a fluctuating semiflexible 
polymer under a compressive force is being carried out by T. Odijk [13]. It is shown 
that thermal fluctuations enhances the bucking threshold in the following way: 



where c is a numerical coefficient of the order one. This can be understood as a 
renormalization of the bending modulus. Thermal fluctuations effectively decrease the 
rigidity of a polymer that is under compressive force. 

To examine the effects of electrostatic interactions in the buckling instability, we 
note that for a charged polymer some new length scales enter into the system. For Debye- 
Huckel electrostatic interaction, two important length scales are the Debye screening 
length and the Bjerrum length Ib = e'^/ekBT, where e measures the charge of an 
electron. We assume that ions with charge Ze are distributed uniformly on the chain 
with minimum distance Zq- Combining these length scales we can achieve the electric 
persistence length for this fluctuating charged polymer that is: Lp^ = {Z"^ / 4){\j:) /Iq^Ib- 
Scaling arguments shows that the total persistence length of a charged chain is given 
by: Lp + Lpe [19j. In a very simplified picture for buckling instability in charged systems 
we can replace the persistence length in the buckling force for a neutral system with the 
total persistence length: 



Based on the dimensional analysis we have considered a general dependence on other 
length scales through a dimensionless function 0. This function captures the physics of 
the buckling instability for a charged system at zero temperature. The above expression 
for the buckling force shows a competition between thermal fluctuations and electrostatic 
repulsion in the buckling instability. Thermal fluctuations decreases the buckling force 
while the repulsion can increase it at flnite temperature. 

The mean flled theory for the buckling instability in charged polymers are 
considered in Ref. [11]. It is shown that the analysis of the buckling problem in 
the charged polymers with nearly rod like geometry are associated with an harmonic 
differential equation with a complicated frequency. Deflning a new function u{s) = 0{s), 
and in the limit of 6{s) <^ 1, the mean fleld equation reads: 




(7) 






(9) 
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Figure 2. Left: Gibbs free energy G, for a buckling polymer versus the perpendicular 
displacement 5 of the polymer at its middle point for different values of the external 
forces /. As one can see, for small forces the system is in mechanical stable equilibrium 
while for large forces the straight state is unstable. Right: the overall stiffness of the 
midpoint is plotted for different values of the applied force. At critical force the stiffness 
vanishes. Simulations are performed for freely rotating ends. 



with 

uj\s) = \j'^ ds'V{\s-s'\). (10) 

The analytical solution of this equation for Debye-Huckel interaction is not possible. 
Here we consider a short range repulsion interaction with a simple form like: 

V{r) = Voe-'/^^. (11) 

In the limit of L ^ Ad, we will have cu^ ~ 2XciVo. Solution to the stability equation, 
reveals that the buckling threshold at zero temperature is: 



fc = n + (12) 

This simply shows that the buckling threshold is linearly increasing with increasing the 
interaction range A^i- 



4. Simulation 



In order to analyze the buckling transition in a charged semifiexible polymer, we simulate 
the system dynamics with ESPResSo software. In descretizing the chain, each individual 
segment has a cylindrical geometry with diameter d = 25nm. For a cylinder with length 
/o and diameter d, the friction coefficients read [2U] : 

ln(/o/rf)' Hk/d)' ^ ' 

In our simulations the time step is dt = 1.2 x 10^^^ s, polymer length L = 10/xm, number 
of monomers = 500 and the viscosity of the solution is ?7 = O.OlPa.s (cytoplasmic 
water). In a typical simulation, the equilibration time was about 3 x 10^ dt. 



Paper submitted to J. Phys. : Condens. Matter 



6 



We first perform tlie simulations for an electrically neutral filament at finite 
temperature. To investigate the stability of a system that is under external compressive 
force /, we study its response to a small transverse displacement of a point in the middle 
of the polymer. To do this favor we apply an external auxiliary force = fj^y to the 
midpoint of the chain. A mechanical work with an amount —f±y{s = L/2) is included to 
the free energy. In this case the total free energy of the system is a function of the applied 
forces: G = G{f,f±). Defining a conjugate variable corresponding to the transverse 
force by, 6 = —dG/df± we can construct the following Legender transformation: 



This is the free energy of the system that is constrained to a constant 5, the midpoint 
displacement. Fig. [2]J^left) shows the numerical results of this free energy for different 
values of the external compressive force /. As one can see for small compressive forces the 
state with 5 = is the stable equilibrium state of the chain. Increasing the compressive 
force and in a critical force fc, the curvature of the free energy changes its sign at 
point 6 = signaling the instability transition. To determine the transition point we 
can expand the free energy for small 6 and define a new variable k, that measures the 
overall stiffness of the midpoint by: 



the transition point can be reached by extrapolating the results to find the state with 
k{fc) = 0. Fig. |2]^right) shows the results for a typical simulation. A numerical 
value of order 1.2pN for the buckling force is achieved for a polymer with stiffness 
Lp = 3mm. This picture for the buckling transition, resembles the physics of the phase 
transition in magnetic systems with changing the external magnetic field. In this case 
it is very natural to ask what happens if one change the temperature, the strength of 
the fluctuations. As in case of magnetic systems we expect a non trivial change in the 
transition point. 

5. Results 

To investigate the effects of thermal fluctuations we first consider an electrically neutral 
filament with bending rigidity that has the same order of magnitude of microtubules. 
For better comparison with cytoplasmic microtubules, we choose a total length of 10 fim 
and apply the pushing force to one of the ends. For boundary conditions we apply the 
clamped condition to both ends. Fig. [3]^left) demonstrates the temperature dependence 
of the critical buckling force for different values of the chain stiffness k. The results at 
very low temperature (T — )• 0), are in agreement with the mean filed result. One can 
see that by increasing the temperature the critical force linearly decreases. 

The temperature dependence of the critical force show different slopes for different 
numerical values of the bending rigidity. To show this behavior in the buckling force 
we investigated the buckling transition at a fixed temperature and vary the rigidity. 



n{f,6)=G-hdG/dU 



(14) 



Q(f,S)=Qoif) + lk{f)6' + 0{6f, 



(15) 
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Figure 3. Left: Critical buckling force versus temperature for a neutral semiflexible 
polymer with two freely rotating ends that is allowed to undergo thermal fluctuations. 
Results are presented for different bending rigidities. Right: Buckling force versus 
bending rigidity is plotted. In both graphs the total length of the polymer is 10 /xm, 
and number of segments in the simulations is n — 500. Numerical error in determining 
the critical force is 0.04 pN . 

Fig. |3]^right) shows the behavior of the critical buckhng force in terms of the bending 
stiffness. This is investigated for different values of the temperature. 

The other issue that we want to address here, is the effects of electrostatic 
interactions in a charged filament immersed in a fluid with mobile charges. For studying 
the electrostatic effects, we use the Deby-Huckel interaction that is already introduced. 
To have a feeling of a real system, we choose the elastic and electrostatic parameters 
of an actin filament with total length L = 0.75/im and persistence length Lp = ISfim. 
Charge density of an actin filament is of order 1 electron in each 2.5A^ and the Bjerrum 
length for water at room temperature is about 0.7nm. Fig H] shows the results of 
the simulations for buckling critical force in terms of the Debye screening length at two 
different temperatures. As one can see, for a range of parameters, there is a competition 
between the electrostatic effects and the effects due to thermal fluctuations. At very 
small Debye screening length (large {Iq/^d)), the critical force reaches the value that 
is already obtained for a neutral filament, increasing the screening length, the critical 
force increases. For a constant screening length thermal fluctuations can either increase 
or decrease the critical force. 

In summary we have considered the effects due to the thermal fluctuations in the 
buckling transition for either a charged or neutral filaments. We have shown that for 
a neutral filament thermal fluctuations can always decrease the buckling critical force 
while for a charged filament the effects of fluctuations can compete with the electrostatic 
effects and either decrease or increase the critical force. 
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Figure 4. Critical buckling force versus Iq/Xd for a semiflexible and charged polymer 
with two freely rotating ends at temperatures 2iK and SOO-ft' . Here Lb ~ 0.7nm, 
K = 7.45 l{)~'^^Nm'^ (for actin filament) and the total length of the polymer is 0.75/U m. 
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